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ABSTRACT 

A  STUDY  OF  BOILING  WATER  PLOW  REGIMES  AT  LOW  PRESSURE 

Mario  P.  Fiori 


Submitted  to  the  Department  of  Nuclear  Engineering  in  partial 

fulfillment  of  the  requirement  for  the  joint  degrees  of 

Nuclear  Engineer  and  Master  of  Science  in  Mechanical  Engineer- 
ing. 

A  comprehensive  experimental  program  to  examine  flow 
regime*  at  pressures  below  100  psia  for  boiling  of  water  in 
tubes  was  carried  out . 

An  electrical  probe,  which  measures  the  resistance  of 
the  fluid  between  the  centerline  of  the  flow  and  the  tube 
wall,  was  used  to  identify  the  various  flow  regimes .  This 
probe  proved  to  be  an  ideal  detection  device,  because  of  its 
simplicity,  reproducibility,  and  accurate  representation  of 
the  flow  pattern  within  the  heated  test  section. 

The  major  flow  regimes  observed  were  bubbly,  slug  and 
annular  flow.  Under  certain  conditions  at  high  flow  rates, 
a  wispy-annular  flow  pattern  was  observed.  The  effects  of 
mass  velocity  (0.2  x  10°  -  2.4  x  1q6  lbm/hr-ft2),  inlet  temp- 
erature (100,  150,  200°F),  exit  pressure  (30,  100  psia), 
quality  (x  =  -10  -  +7  percent),   purity  (9,  40  PPM  NaCl;  1-3 
megohm-cm),  length  (L/D=30,  60,  90),  diameter  ( 0.094,  0.242 
in.),  and  orientation  (vertical  and  horizontal)  on  the  flow 
regimes  were  studied.  Flow  regime  maps  on  coordinates  of 
mass  velocity  and  quality  are  presented  for  these  conditions. 

Bubbly  and  slug  flow  occurred  primarily  in  the  sub- 
cooled  region ,  while  fully  developed  annular  flow  was  reached 
at  equilibrium  qualities  between  2  and  4  percent 0   The  tran- 
sitions between  the  different  flows  were  shifted  to  regions 
of  increased  subcooling  when  velocity,  pressure,  and  heat 
flux  increased,  and  when  inlet  temperature  decreased.   Purity 
and  geometry  had  little  affect  on  the  flow  regime  boundaries. 

The  shifting  of  the  transitions  is  related  to  the  agglom- 
eration point,  which  Is  that  point  at  which  the  bubbles  so 
coalesce  that  slug  flow  is  first  observed.  The  agglomeration 
point  depends  on  the  point  of  incipient  boiling,  the  number 
of  bubbles  in  the  flow,  and  the  number  of  collisions  per 
bubble.  These  latter  quantities  in  turn  depend  on  velocity, 
temperature,  pressure,  and  heat  flux. 

The  flow  regime  information  obtained  in  this  study  should 
be  of  value  in  correlating  and  interpreting  low  pressure  heat- 
transfer  data.  The  flow  regime  data  were  found  to  be  useful 
in  explaining  the  effect  of  inlet  temperature  on  burnout  heat 
flux. 

Thesis  Supervisor:   Arthur  E.  Bergles 

Title:  Assistant  Professor  of  Mechanical  Engineering 
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NOMENCLATURE 
A      =  heat  transfer  area  (ft2) 


c      =  specific  heat  (Btu/lbm°P) 


D  =  diameter  (in.) 

G  =  mass  velocity  (lbra/hr-ft2) 

h  =  heat  transfer  coefficient  (Btu/hr-ft2oF) 

he  =  enthalpy  at  exit  (Btu/lbm) 

hf  =  enthalpy  of  saturated  liquid  (Btu/lbm) 

hf  =  latent  heat  (Btu/lbm) 

h.  =  enthalpy  at  entrance  to  test  section  (Btu/lbm) 

I  =  shunt  voltage  (mv) 

k  =  thermal  conductivity  (Btu/hr-ft °F) 

L  =  length  (in.,  ft) 

F  =  pressure  (psia) 

P  =  pressure  at  exit  of  test  section  (psia) 

q  =  rate  of  heat  flow  (Btu/hr) 

q/A  =  heat  flux  (Btu/hr -ft2) 

(q/A).  =  heat  flux  at  incipient  boiling  (Btu/hr-ft2) 

(q/A)cr=  critical  heat  flux  (Btu/hr-ft2) 

Rj  =  oscilloscope  internal  impedance  (ohms) 

R  =  resistance  between  probe  and  tube  wall  (ohms) 

R  =  variable  resistance  (ohms) 

T.  =  water  temperature  at  inlet  or  exit  of  tube  (°F) 


Th     =  water  bulk  temperature  (°F) 


b 
T      =  wall  temperature  (°F) 

T      =  saturation  temperature  (°F) 
s 


=  voltage  across  test  section  (volts) 
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VR  =  battery  voltage  (volts) 
VT  =  voltage  across  variable  resistor  (volts) 
w  =  mass  flow  rate  (lbm/hr) 
x   =  quality 
[i     =  viscosity  (lbm/hr-ft) 

Flow  Regime  Designations 

B   =  #  =  bubbly 

BTS  «  #  =  bubbly-to-slug  transition 

S   =  O  =  slug 

STA  =  [J  =  first  indication  of  slug-to-annular  transition 

STA  =  A  =  slug -to-annular  transition 

A   =  ■  =  annular 

WA      =  Jfy  =  wispy  annular 
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Chapter  I 
INTRODUCTION 
1 .1  Background  of  Problem 

In  recent  years  the  interest  in  two-phase,  gas-liquid 
flow  has  dramatically  expanded  as  is  indicated  by  the  tre- 
mendous increase  in  the  number  of  experiments  and  published 
papers  in  this  field.   If  one  examines  the  literature,  one 
will  find  that  a  relatively  large  number  of  models  have 
been  proposed  for  analyzing  and  correlating  experimental 
pressure-drop,  heat-transfer,  critical  heat  flux,  and  void- 
fraction  data.  These  models  frequently  depend  on  the  flow 
regimes  (or  flow  patterns)  which  are  geometric  descriptions 

of  the  distribution  of  the  liquid  and  gas  phases  within  a 

( 1  )* 
channel.  Gouse*  '   considers  knowledge  of  the  flow  regime 

to  be  as  important  as  knowing  whether  the  flow  is  laminar 

or  turbulent  in  single-phase  flow. 

( 2) 

Bakerv  '  showed  that  it  was  not  practical  to  disregard 

the  flow  regimes  in  devising  a  two-phase  pressure-drop 

( 3) 

correlation.  Martinelli's^  '  assumptions  in  his  pressure- 
drop  correlations  tend  to  limit  them  to  annular  flow. 
Lavin  and  Young'  '  propose  boiling  heat-transfer  correla- 
tions for  refrigerants,  which  are  dependent  on  the  flow 
regimes.  Numerous  models  for  the  critical  heat  flux  condi- 
tion are  based  on  the  postulate  of  spray  annular  flow; 
however,  there  are  some  analyses  which  assume  that  fog  flow 
is  present . 

♦Numbers  in  parenthesis  refer  to  References  on  p.  80 . 
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It  is  apparent,  then,  that  it  is  important  to  know 
the  flow  regime  if.  one  is  to  utilize  an  appropriate  model; 
however,  direct  relationships  between  flow  regimes  and 
observables,  i.e.,  pressure  drop,  heat-transfer  coeffi- 
cient, critical  heat  flux,  and  void  fraction  have  not  been 
clearly  established.^   it » still  appears,  though,  that 
flow  regime  observations  will  be  helpful  in  reducing  the 
complications  of  correlating  the  quantities  of  basic 
Interest . 

A  survey  of  the  literature  indicates  that  -flow  regime 
studies  have  almost  always  been  conducted  under  adlabatic 
conditions.  However,  exploratory  studies  have  shown  that 
flow  patterns  are  considerably  altered  when  there  is  heat 
addition,  particularly  in  low  pressure  systems.^  ' 
Comprehensive  flow  regime  maps  have  been  prepared  at 
Dynatech^''*  '  and  Harwell^-)  for  high  pressure  boiling 
water  systems.  No  similar  study  has  been  conducted  at  low 
pressures  (below  100  psia). 

Low  pressure  systems  are  receiving  increased  atten- 
tion, particularly  in  the  nuclear  reactor  and  desaliniza- 
tion  programs.  In  a  recent  study  of  critical  heat  fluxes 
at  low  pressure,  Lopina^   '  observed  a  pronounced  effect 
of  inlet  temperature.  For  constant  geometry,  flow  rate, 
and  pressure  level,  the  critical  heat  flux  Increased  as 
the  inlet  temperature  increased.  This  effect,  which  is 
contrary  to  what  is  observed  at  high  pressures,  appears  to 
be  related  to  the  flow  regimes.  The  phenomenon  of 
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upstream  boiling  burnout  in  uniformly  heated  channels  also 
appears  to  be  related  to  flow  regimes,  in  particular,  the 
slug  flow  regime  . 

1 .2  Flow  Regimes 

There  have  been  various  descriptions  of  flow  regimes 
by  different  investigators .  These  classifications  have 
ofter  been  quite  arbitrary,  and  are  also  complicated  by  the 

fact  that  vertical  and  horizontal  flows  cause  different 

(2) 

regimes.  Bakerv  ;  described  no  less  than  seven  different 

regimes  in  adiabatic  two-phase,  gas-liquid  flow  in  horizon- 
tal pipes.  Vohr'   '  and  Kepple  and  Tung'   ^  have  presented 
complete  bibliographies  on  studies  of  adiabatic  flow 
patterns . 

The  flow  regimes  which  may  be  encountered  in  flows  in 
tubes  with  vertical  upflow  can  be  described  in  the  order  of 
increasing  gas  or  vapor  content.   If  very  little  gas  or 
vapor  is  involved,  the  flow  will  be  characterized  by  isola- 
ted vapor  or  gas  bubbles  dispersed  rather  uniformly  in  the 
continuous  liquid  phase,  This  is  known  as  "bubbly  flow". 
More  bubbles  appear  with  an  increase  of  vapor  or  gas  until 
the  bubbles  begin  to  touch  one  another  and  agglomerate. 
As  more  bubbles  agglomerate,  the  gas  begins  to  flow  in 
long  cylindrical,  bullet-shaped  bubbles,  Between  the 
bubbles  are  slugs  of  liquid  and  hence  "slug  flow". 
Further  increase  of  the  gas  content  destroys  the  stability 
of  the  bubbles,  producing  "churn  flow"  or  "semi  -annular 
flow"  in  which  slugs  of  liquid  separate  adjacent  gas 
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bubbles  only  temporarily.  As  the  gas  content  continues  to 
increase,  the  liquid  flows  on  the  wall  as  an  annular  film 
and  the  gas  or  vapor  flows  in  the  center  of  the  tube.  This 
geometry  is  called  "annular  flow".   "Fog  flow"  or  "dispersed 
flow"  then  follows  where  the  gas  or  vapor  phase  occupies 
the  entire  tube  but  small  drops  of  liquid  are  dispersed 
within  it.  This  flow  is  also  called  "mist  flow". 

Other  flow  patterns  which  have  been  encountered  are 
"froth  flow",  "spray  annular  flow"  and  "wispy  annular  flow". 
Froth  flow  is  quite  similar  to  bubbly  flow  except  the 
bubbles  are  so  numerous  that  the  flow  has  a  frothy  appear- 
ance. Spray-annular  flow  is  a  combination  of  mist  and  annu- 
lar flow.  Some  of  the  liquid  flows  in  an  annular  film  and 
some  flows  as  small  droplets  in  the  vapor  core.  Wispy  annu- 
lar is  similar  to  spray  annular  except  the  flow  in  the  core 
consists  of  agglomerated  liquid  instead  of  small  droplets. 

There  are  very  similar  flow  descriptions  for  horizontal 
flow.  However,  stratification  is  sometimes  encountered  in 
horizontal  flow  and  several  new  regimes  are  introduced. 
"Plug  flow"  is  similar  to  slug  flow  except  that  the  bubbles 
move  along  the  upper  part  of  the  channel  or  pipe .   In 
"stratified  flow"  the  liquid  flows  at  the  bottom,  and  gas  or 
vapor  at  the  top  of  the  pipe. 

These  regimes  appear  to  be  quite  well  defined.  However, 
for  practical  reasons  it  is  desirable  to  minimize  the 
number  of  regimes.  Gouse*  '  suggests  that  there  are  only 
four  distinct  flow  patterns  for  vertical  upflow:   bubbly, 
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slug,  annular,  and  mist;  while  all  the  others  which  have 
been  observed  are  merely  transitions  from  one  to  another. 

The  above  flows  were  described  for  adiabatic  condi- 
tions. In  heated  systems  such  clean  geometric  descriptions 
of  the  flew  are  not  expected  because  of  several  character- 
istics not  encountered  under  adiabatic  conditions.  Some  of 
these  characteristics  are: 

1)  Vapor  is  continually  formed  along  the  tube.  This 
creates  a  volume  flow  rate  which  continuously  increases 
along  the  tube;  thus  the  flow  regimes  would  be  expected  to 
be  different  in  different  sections  of  the tube . 

2)  Bubbles  may  be  generated  and  disturb  the  flow 
pattern. 

In  any  case,  the  adiabatic  flow  regimes  are  a  basis  for 
interpreting  and  classifying  the  flow  regimes  with  heat 
addition.  The  classification  and  both  qualitative  and  semi- 
quantitative description  of  the  flow  patterns  observed  in 
this  study  are  discussed  in  Section  2.3.3. 

1.3  Flow  Regime  Detection 

Numerous  methods  of  studying  flow  regimes  in  heated 
systems  have  been  developed.  High-speed  still  and  movie 
photography  has  been  most  popular.  More  exotic  methods 
which  have  been  used  include  isokinetic  flow  sampling, 
gamma  or  beta  ray  attenuation,  and  hot-wire  anemomet-ry. 
Each  of  these  methods  has  drawbacks.  Either  there  is  too 
much  data  to  analyze,  as  in  the  case  of  movies,  or  the 
methods  are  Very  complicated  and  expensive .  The  electri- 
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cal  resistance  probe  which  was  developed  by  Solomon^13' 
specifically  for  flow  regime  detection,  is  an  inexpensive, 
simple,  and  reliable  device. 

This  probe  measures  the  resistance  between  the  center- 
line  of  the  fluid  flow  and  the  tube  wall.  A  bridge  of 
liquid  across  the  tube  will  cause  a  closed  circuit,  and 
gas  or  steam  across  the  tube  will  result  in  an  open  circuit 
Nassos^   '  and  Haberstroh^  ^'   have  successfully  used  this 
type  of  probe  for  adiabatic  flow.  Griffith'1  '   adapted 
this  probe  for  a  high  pressure  boiling  water  system  and 
delineated  the  slug- to -annular  transition.  Dynatech' ' ''  ' 
has  recently  been  using  this  technique  to  map  boiling  water 
flow  regimes  at  high  pressures  ( 500-1000  psia). 

1 .4  Scope  of  Research 

A  comprehensive  study  to  examine  boiling  water  flow 
regimes  in  tubes  of  small  diameter  at  pressures  below  100 
psia  was  undertaken.  A  resistance  probe  described  in 
detail  in  Section  2.2  was  chosen  as  the  primary  means  of 
flow  regime  detection;  however,  a  visual  section  was  also 
employed.  The  experimental  program  consisted  of  an  inves- 
tigation of  the  effects  of  velocity,  inlet  temperature, 
exit  pressure,  quality,  (heat  flux),  purity,  length,  diam- 
eter, and  orientation  on  the  flow  regimes.  The  first 
phase  of  the  study  was  devoted  to  probe  development  and 
interpretation. 
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Chapter  II 
EXPERIMENTAL  PROGRAM 
2.1  Description  of  Apparatus 

2.1.1  Hydraulic  System 

The  experimental  facility  used  was  the  low-pressure 
test  loop  located  in  the  MIT  Heat  Transfer  Laboratory.  The 
basic  apparatus  was  designed  and  constructed  in  196l.^1^' 

A  schematic  of  the  loop  is  presented  in  Fig.  1.  The 
pipings  and  fittings  are  made  of  brass  and  stainless  steel 
for  corrosion  resistance.  Rayon  reinforced  rubber  hose  is 
used  where  flexible  connections  are  required.  Distilled 
water  is  circulated  by  a  bronze,  two-stage,  regenerative 
pump  providing  a  discharge  pressure  of  260  psig  at  3.6  gpm. 
The  pump  is  driven  through  a  flexible  coupling  by  a  3-hp 
Allis-Chalmers  induction  motor.  A  Pulflo  filter  is  in- 
stalled at  the  pump  inlet.   Pressure  fluctuations  at  the 
outlet  of  the  pump  are  damped  out  by  means  of  a  2.5-gal 
Greer  accumulator  charged  with  nitrogen  to  an  initial 
pressure  of  40  psig.  This  accumulator  contains  a  flexible 
bladder-type  separator  which  prevents  the  nitrogen  from 
being  absorbed  by  the  system  water.  After  the  accumulator, 
the  flow  splits  into  the  by-pass  line  and  the  test-section 
line. 

In  the  test-section  line,  fluid  flows  through  a 
Fischer-Porter  flowrator  followed  by  a  preheater,  thence 
through  a  Hoke  metering  valve  and  the  test  section,  after 
which  it  merges  with  fluid  from  the  by-pass  line.  The 
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flow  then  goes  through  the  heat  exchanger  and  returns  to 
the  pump.  The  preheater  consists  of  four  Chromalox  heaters 
of  approximately  6  kw  each.  Three  of  these  are  controlled 
simply  with  "off-on"  switches  while  the  fourth  can  provide 
a  continuous  range  from  0  to  6  kw  by  means  of  a  bank  of 
two  variacs  mounted  on  the  test  bench.  Quick-action  James- 
bury  ball  valves  are  installed  before  the  inlet  to  the 
flowrator  and  after  the  exit  from  the  test  section.  This 
permits  quick  isolation  of  the  test  section  to  minimize 
fluid  loss  when  conducting  burnout  tests.  The  exit  valve 
is  also  used  to  adjust  the  test-section  pressure. 

Flow  through  the  by-pass  line  is  controlled  by  a  ball 
valve  on  each  side  of  which  there  is  a  300-psig  pressure 
gage.   Pump  operating  pressure,  and  hence  the  pressure  up- 
stream of  the  test  section,  is  controlled  by  this  valve . 

The  heat  exchanger  is  a  counterflow  type  with  system 
water  flowing  in  the  inner  tube  and  city  water  in  the 
outer  annulus .  Except  at  very  high  power  levels,  the  heat 
exchanger  maintained  a  constant  pump  inlet  temperature.  A 
Fulflo  filter  is  installed  on  the  city  water  line  to 
reduce  scale  formation  in  the  exchanger. 

The  distilled  water  was  de-ionized  continuously, 
except  during  contamination  studies,  by  passing  a  portion 
of  the  flow  through  four  mixed-bed  resin  demineralizer 
cartridges  installed  in  a  Barnstead-  "Bantam"  demineralizer 
unit.  A  4.7-gal  degassing  tank  is  provided  with  five  elec- 
trical heaters  (3-200  vac  and  2-110  vac).  This  tank  also 
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serves  as  a  surge  tank.  A  15-gal  stainless-steel  storage 
tank  for  filling  the  system  is  mounted  directly  above  the 
degassing  tank  and  can  be  filled  with  distilled  water  from 
standard  5-gal  bottles  with  a  small  Hypro  pump. 

2.1.2  Power  Supply 

Test-section  electrical  power  is  provided  by  two 
36-kv;  dc  generators  connected  in  series.  Each  generator  is 
rated  at  12  volts  and  3000  amperes.  The  power  control  con- 
sole permits  coarse  or  fine  control  from  0  to  2b   volts.  A 
water-cooled  shunt  installed  in  parallel  with  the  test  sec- 
tion protects  the  generators  against  the  shock  of  the  sudden 
open  circuit  which  occurs  at  burnout.   Power  is  transmitted 
from  the  main  bus  to  the  test  section  by  water-cooled  power 
leads.   In  the  present  tests,  the  downstream  or  exit  power 
lead  was  at  ground  potential.  Rubber  hose  connected  both 
inlet  and  exit  chamber  plenums  to  the  main  loop  in  order 
to  electrically  isolate  the  test  section.  The  exit  of  the 
test  section  was  then  separately  grounded  (see  schematic 
of  test  section,  Pig.  2). 

2.1.3  Inst  rume  nt  a  t  i  on 

All  temperatures  were  measured  by  copper-constantan 
thermocouples  made  from  30-gauge  Leeds  and  Northrup  duplex 
wire.  The  test-section  inlet  temperature  was  measured  by 
a  thermocouple  directly  in  the  fluid  stream,  upstream  of 
the  Hoke  metering  valve.  The  thermocouple  was  introduced 
at  this  location  through  a  Conax  fitting  equipped  with  a 
lava  sealant. 
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All  pressures  were  read  on  Bourdon-type  gages  located 
as  shown  in  Fig.  1.  The  test-section  inlet  and  exit  pressures 
were  measured  with  Helicoid  8-1/2  in.  gages  of  200  psig  and 
100  psig,  respectively.  Both  are  specified  to  an  accuracy 
of  +0.25$  of  full  scale. 

A  variety  of  metering  tubes  and  floats,  which  could  be 
installed  interchangeably  in  the  basic  Fischer-Porter  flow- 
meter housing,  provided  measurement  of  the  test-section  flow 
from  1.5  to  4000  lbm/hr . 

The  voltage  drop  across  the  test  section  was  read 
directly  on  a  Weston  multiple -range  dc  voltmeter  with  a  spec- 
ified accuracy  of  +  1/2  percent,   The  current  flow  was 
determined  by  using  the  Minneapolis-Honeywell,  Brown  recor- 
der to  measure  the  voltage  drop  across  a  calibrated  shunt 
(60.17  amp/mv)  in  series  with  the  test  section. 

The  details  of  the  resistance  probe  are  discussed  in 
Section  2.2.   Polaroid  cameras  were  used  for  the  photo- 
graphic study.  The  camera  settings  and  procedures  are  dis- 
cussed in  Appendix  B. 

2.1.4  Operating  Procedures 

After  the  test  section  was  installed,  the  loop  and  de- 
gassing tank  were  filled  with  distilled  water  from  the 
supply  tank.  The  water  in  the  degassing  tank  was  then 
brought  to  a  boil  while  the  loop  water  was  circulated  with 
the  heat-exchanger  coolant  off.  Degassing  was  accomplished 
by  by-passing  a  portion  of  the  loop  water  into  the  top  of 
the  vigorously  boiling  degassing  tank.  This  was  continued 
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until  the  temperature  of  the  loop  rose  to  approximately 
180°F.  A  standard  Winkler  analysis  described  in  Ref .  17 
indicated  that  this  method  of  degassing  reduced  the  air 
content  to  less  than  0.1  cc  air/liter.  After  completing 
degassing,  the  loop  water  was  adjusted  to  the  desired  in- 
let temperature  and  the  desired  pressures  were  set.  A 
pressure  of  from  210  psig  to  260  psig  was  maintained  at  all 
times  on  the  upstream  side  of  the  Hoke  metering  valve  to 
prevent  system-induced  instabilities.  The  generators  were 
then  started  and  allowed  to  warm  up. 

Nine  burnout  tests  were  made  in  order  to  have  reliable 
critical  heat  flux  data  for  conditions  similar  to  the  planned 
flow  regime  studies.  This  knowledge  aided  in  avoiding  burn- 
out of  the  probe  test  sections.  Power  was  applied  to  the  test 
section  in  small  increments  while  maintaining  a  constant  flow 
rate  and  inlet  temperature.  The  test-section  flow  rate, 
shunt  voltage  (test-section  current),  inlet  temperature,  and 
exit  pressure  were  recorded  for  each  voltage  setting.  The 
critical  or  burnout  heat  flux  was  defined  to  be  that  heat  flux 
which  caused  physical  failure  of  the  test  section.  After 
commencing  flow  regime  studies,  it  was  seen  that  the  area  of 
interest  was  usually  far  removed  from  (q/A)cr.  Consequently 
no  further  burnout  studies  were  made.  Appendix  A  contains  a 
tabulation  of  the  burnout  points. 

The  same  data  were  recorded  for  the  flow  regime  runs. 
A  run  consisted  of  taking  data  at  a  particular  mass  flow 
rate  for  a  given  set  of  conditions,  i.e.  D,  L/D,  T±,    Fq . 
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As  the  power  was  increased,  various  responses  on  the  oscil- 
loscope indicated  the  type  of  flow  (see  Probe  Interpretation, 
Section  2.3).  Water  purity  was  checked  by  means  of  a 
Barnstead  Purity  Meter  (0.1  -  18  megohm-cm)  or  an  Industrial 
Instruments  Solu-Bridge  Meter  (0  -  30  PPM  NaCl). 

2.2  Resistance  Probe  and  Test  Section 

2.2.1  Principle  of  Operation 

The  resistance  probe  was  used  to  determine  the  flow 
regime  at  the  exit  of  the  test  section.  The  probe  is  a  wire 
placed  in  the  center  of  the  tube  as  shown  in  Fig.  2.  This 
wire  is  covered  by  teflon  sleeving  except  for  the  tip  of  the 
wire.  The  circuit  diagram  for  the  probe  is  shown  in  Fig.  3a. 
When  there  is  a  bridge  of  water  across  the  tube  at  the  test- 
section  exit,  there  is  essentially  a  closed  circuit  between 
the  probe  and  the  tube  wall.  This  results  in  a  maximum  volt- 
age across  resistor  R,  which  is  measured  by  an  oscillo- 
scope. On  the  other  hand,  steam  passing  between  the  probe 
and  the  wall  resembles  an  open  circuit  giving  a  voltage 
reading  across  FL.  near  zero. 

2.2.2  Probe  Integration  with  Test  Section 

Two  types  of  test  sections  were  used  for  this  study. 
The  first  type  was  used  for  preliminary  burnout  studies. 
The  second  type  incorporated  the  resistance  probe  and  a 
sight  section.  The  sight  section  was  not  always  used 
because  the  rubber  gasket  would  not  hold  for  exit  pressures 
above  40  psig.  The  heated  sections  were  made  of  standard 
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304  stainless  tubes.  The  burnout  test  sections  and  the  first 
probe  test  section  had  an  0.242-in.  i.d.  and  0.3125-ln.  o.d. 
while  the  second  probe  test  section  had  an  i.d.  of  0.094  in. 
and  o.d.  of  0.120  in.  The  details  of  the  two  test  sections 
are  shown  in  Fig.  2. 

In  the  burnout  test  section,  3/4- in.  brass  bushings 
were  cilver-soldered  to  the  stainless  steel  tube  to  form 
the  power  connections  upstream  and  downstream.  A  pressure 
tap  was  installed  in  each  exit  bushing  with  a  No.  80, 
0.0135-in.  diameter  hole  drilled  through  the  tube  wall. 
The  pressure  tap  was  placed  here  rather  than  in  the  exit 
plenum  chamber  because  of  the  choking  phenomena  whereby  the 
plenum  pressure  is  lower  than  the  exit  pressure.^   ' 

The  probe  test  section  had  pressure  taps  located  at 
both  the  inlet  and  exit.  The  exit  section  was  modified  to 
include  the  probe  and  pyrex  sight  section.  The  Inside 
diameter  of  the  glass  tubing  was  equal  to  the  outside  diam- 
eter of  the  test  section.  This  caused  a  slight  expansion 
at  the  entrance  to  the  visual  section. 

The  brass  supporting  structure  for  the  probe,  pressure 
tap,  and  sight  section  also  served  as  the  exit  power  connec- 
tion. The  pressure  tap  is  similar  to  the  previous  one 
described. 

The  probe  was  fabricated  from  0.035-in.  diameter  temp- 
ered stainless-steel  spring  wire  for  the  0.242-in.  diameter 
tube  and  from  0.015-in.  wire  for  the  0.094-in.  diameter 
tube.  Teflon  spaghetti  was  slipped  over  this  wire  and 
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stretched  until  the  teflon  set  snugly  on  the  wire  .   The 
final  probe  diameters  were  0.060  in.  and  0.031  in.  Even  the 
small  diameter  probe  had  adequate  stiffness  to  prevent  vi- 
bration which  could  create  spurious  signals. 

The  straight  probe  was  pushed  into  the  test  section 
from  the  packing  gland  side  until  it  touched  the  opposite 
wall.  The  tip  was  then  carefully  bent  and  the  teflon  was 
forced  around  the  wire  until  only  a  1/4  in.  of  probe  tip  was 
bare . 

The  teflon  acts  as  an  insulation  between  the  probe  and 
the  tube  wall.  More  important,  since  teflon  is  not  wetted 
by  water,  there  is  no  chance  of  a  short-circuit  between  the 
exposed  end  of  the  wire  and  the  tube  wall.  The  drying  time 


of  the  teflon  after  passage  of  a  slug  of  liquid  is  extremely 
les: 

(15) 


_3 
short,  less  than  10   seconds  as  observed  with  an  oscillo- 


scope 

2.2.3   Probe  Sensitivity 

A  schematic  diagram  of  the  probe  circuit  is  shown  in 
Fig.  3a.   It  can  be  shown  that  the  voltage  reading,  V-j-,  is 
related  to  the  resistance  of  the  water-steam  mixture  in  the 
test  section,  R  ,  by  the  equation 

v.-   ,t  Vfi1 (D 

1  Rw(^ +  £>  + 1 

Furthermore,  it  can  be  shown  that  the  maximum  sensitivity, 
defined  by  a  maximum  in  (3^)  is  obtained  when 

i-i+A  (2) 
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At  this  maximum  sensitivity,  the  instrument  voltage,  V,.,  is 

related  to  the  battery  voltage,  V^,  by  the  equation 

.b 

VR 

?!--?-  (3) 

The  dependence  of  VT  on  R  ,  for  VV,  =  6  volts  and  RT  = 

10  ohms,  is  indicated  in  Fig.  3b.   The  lowest  value  of  R 

occurs  when  saturated  liquid  is  flowing  in  the  test  section. 

During  the  course  of  the  experiments,,  R  for  saturated  liq- 

w 
r  6 

uid  was  found  to  vary  between  10^  and  2  x  10  ohms,  For  the 

lower  limit  of  R  =  10^  ohms,  it  can  be  seen  that  for  R  = 

10^  ohms,  the  maximum  sensitivity  will  occur  in  the  liquid 

region.  For  R  =  oo  the  maximum  sensitivity  will  occur  at  a 

much  higher  R  corresponding  to  a  flow  pattern  with  a  higher 

void  fraction .  Consequently,  by  adjusting  R  ,  the  maximum 

sensitivity  can  be  chosen  to  correspond  to  a  flow  regime 

boundary.  However,  since  these  boundaries  are  not  known  in 

advance,  it  is  desirable  to  maintain  R  at  a  fixed  value  for 

each  run.   It  wa3  found  most  convenient  to  adjust  R  so  that 

maximum  sensitivity  was  obtained  in  the  bubbly  region.  This 

could  not  always  be  achieved  due  to  the  restriction  imposed 

by  the  internal  resistance  of  the  oscilloscope  .  This,  however, 

did  not  affect  the  results  to  any  noticeable  degree  since  the 

variable  of  interest  was  the  shape  rather  than  the  value  of  Vj. 

2.3  Probe  Interpretation 

The  output  of  the  probe  for  typical  runs  is  described 
below.  Throughout  the  interpretation  phase  of  the  experi- 
ment, Polaroid  photos  were  simultaneously  taken  of  the  exit 
flow  and  of  the  oscilloscope  trace. 
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2.3.1  Velocity  Effect  on  Probe  Signal  Amplitude 
With  stationary  water  in  the  test  section,  a  large 

signal  was  obtained.  When  the  pump  was  started  and  water 
began  to  flow,  the  signal  decreased  by  one  to  two  volts. 
As  the  water  velocity  increased,  the  voltage  decreased, 
although  the  decrease  was  much  less  than  that  observed  at 
startup.   Other  experimenters^  -?'  also  observed  this  effect 
and  speculated  that  this  is  an  "electroklnetic"  effect, 
associated  with  the  charge-carrying  capacity  of  the  water. 

2.3.2  Temperature  Effect  on  Probe  Signal  Amplitude 
When  the  power  was  applied  to  the  test  section,  the 

bulk  fluid  temperature  increased  and  the  voltage  trace  in- 
creased. This  voltage  increased  to  a  maximum  as  the  satu- 
ration temperature  was  reached.  This  was  expected  since 
the  conductivity  of  the  water  increases  with  temperature, 
causing  a  greater  voltage  drop  across  the  variable  resis- 
tor, R  .  To  Insure  that  this  effect  was  only  due  to  temp- 
erature, power  was  shut  off  and  the  water  temperature  was 
increased  by  means  of  the  preheaters0  The  same  effect  was 
observed. 

2.3.3  Flow  Regimes  Related  to  Probe  Signal 

The  number  of  distinct  flow  regimes  of  two  phase  flow 
was  shown  to  be  somewhat  arbitrary.  Four  main  flow 
patterns  were  recognized  in  the  present  experiments: 
bubbly,  slug,  annular,  and  wispy  annular  flow.  The  physi- 
cal or  geometric  description  of  the  observed  flow  patterns 
is  very  similar  to  that  given  in  Section  1.2. 
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Bubbly  flow  is  characterized  by  a  continuous  liquid 
phase  with  small  vapor  bubbles  dispersed  somewhat  uniformly 
throughout  the  tube.   The  characteristic  bullet-shaped 
bubbles  of  slug  flow  could  not  be  identified  under  the 
present  high  heat  flux  conditions .  Slug  flow  was  recog- 
nized whenever  large,  irregular  bubbles,  whose  size  was  on 
the  order  of  the  diameter  of  the  tube  had  somewhat  continu- 
ous slugs  of  water  between  them.   Annular  flow  consisted  of 
a  thin  film  of  liquid,  usually  quite  turbulent,  flowing  along 
the  wall,  while  vapor  and  fine  droplets  flowed  in  the  center 
of  the  tube.  This  has  often  been  called  spray-annular  flow, 
but  for  the  purposes  of  this  report,  annular  will  be  used. 
The  character  of  the  core  flow  was  Inferred  from  the  charge 
carrying  spray  effect  which  is  discussed  in  Sec.  2.3.4. 
Wispy  annular  flow  was  characterized  by  a  continuously  in- 
creasing void  fraction  without  the  distinct  characteristics 
of  the  above  three  flows. 

Since  the  classification  of  these  flow  regimes  is 
rather  qualitative,  different  individuals  may  interpret 
the  flows  differently.  The  data  in  Figs .  4,  5,  and  6  were 
chosen  so  as  to  indicate  the  semi-quantitative  definitions 
of  the  flows  which  were  used  throughout  this  study. 
Figure  7  presents  a  summary  of  the  probe  traces  for  the 
various  flow  patterns. 

Bubbly  Flow  In  Photo  I  in  Fig,  4,  a  2.8-volt  signal  with 
very  small  fluctuations  is  observed .  These  small  fluctua- 
tions are  due  to  bubbles  which  can  be  seen  in  the  sight 
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section.  Figure  7a  shows  the  schematic  of  the  oscillo- 
scope trace  for  bubbly  flow.   Incipient  boiling  could  not 
be  detected  by  the  probe  due  to  the  small  amount  of  vapor 
generation.   It  could  not  be  detected  with  the  sight  sec- 
tion either,  because  these  bubbles  recondensed  in  the 
subcooled  bulk  flow  before  visual  observation  could  be 
made  . 

Bubbly-to-Slug  Transition  The  transition  from  bubbly  to 
slug  (BTS)  flow  was  first  defined  when  the  small  pips 
previously  seen  in  bubbly  flow  became  large  enough  to  pro- 
duce a  1 .0  -  1.5-volt  deflection.  Photo  II  graphically 
shows  a  large  bubble  just  leaving  the  heated  section. 
This  bubble  caused  a  1.3-volt  deflection  on  the  scope 
trace.  This  transition, schematically  shown  in  Fig.  7b, 
was  well  defined  and  reproducible . 

Slug  Flow  Slug  flow  was  considered  to  be  present  when 
the  probe  signal  alternated  between  all-liquid  and  all- 
vapor  levels.   In  Photo  III  it  is  noted  that  clean 
bubbles  and  slugs,  of  the  variety  found  in  adiabatic  flow, 
do  not  exist .  The  bubble  is  like  a  vapor  void  without  a 
definite  head  , or  tall,  and  the  slugs  have  small  bubbles 
dispersed  throughout.  At  high  flow  rates,  the  slug  trace 
was  somewhat  different,  as  seen  in  Photo  IV.  The  frequency 
of  the  full  deflection  is  increased  dramatically  and  the 
flow  appears  to  be  very  agitated.  Figures  7c  and  7d  show 
the  first  full  slug  and  slug  flow,  respectively. 
Slug-to -Annular  Transition  The  transition  criterion  for 
the  slug  to  annular  (STA)  was  more  subjective  than  the 
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other  criteria,,   From  Fig.  7d  it  is  seen  that  for  slug  flow, 
the  voltage  deflections  average  out  between  (VT)    and  0. 
In  Fig.  7e  it  appears  that  the  slugs  are  thinner  and  a 
larger  percentage  of  vapor  is  passing  the  probe.  The  STA 
transition  was  defined  when  the  voltage  deflections  appeared 
to  rise  from  0  volts  up,  rather  than  from  (V-j.)    down.   An 
example  of  this  flow  is  seen  in  Photo  V0  It  is  noted  that 
it  would  be  very  difficult  to  distinguish  this  flow  from 
annular  flow,  particularly  at  high  G,  using  the  sight  sec- 
tion alone.  As  test-section  power  increases  further 
(Photo  VI,  Fig.  If),    the  voltage  trace  is  similar  to  that  of 
the  initial  description  of  the  STA  transition  except  that 
the  peak  voltages  are  not  full  deflections.  This  flow  is 
still  the  transitional  flow  from  slug  to  annular.   Physi- 
cally, it  appears  as  if  most  of  the  tube  consists  of  annular 
flow,  but  much  of  the  liquid  transport  is  still  in  the  very 
rapidly  moving  slugs . 

Annular  Flow  With  greater  power,  the  voltage  deflections 
become  smaller,  thus  indicating  that  the  annular  flow  region 
is  being  approached.  Annular  flow  was  defined  to  be  that 
flow  which  would  give  an  average  deflection  no  greater  than 
0.5  volts  (Fig.  7g).  This  definition  did  allow  larger  de- 
flections to  occur,  but  the  frequency  of  these  large  deflec- 
tions was  greatly  reduced.   At  low  velocity  annular  flow 
was  well  defined  since  the  voltage  was  practically  zero 
(Photo  VII).   Photo  VIII  is  another  example  of  annular  flow. 
The  increased  velocity  of  the  vapor  core  caused  large  wave- 
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other  criteria.   From  Fig.  7d  it  is  seen  that  for  slug  flow, 

the  voltage  deflections  average  out  between  (Vj)  a  and  0. 

In  Fig.  7e   it  appears  that  the  slugs  are  thinner  and  a 

larger  percentage  of  vapor  is  passing  the  probe.  The  STA 

transition  was  defined  when  the  voltage  deflections  appeared 

to  rise  from  0  volts  up, rather  than  from  (Vx)  „   down.   An 

3  v  I 'max 

example  of  this  flow  is  seen  in  Photo  V.   It  is  noted  that 
it  would  be  very  difficult  to  distinguish  this  flow  from 
annular  flow,  particularly  at  high  G,  using  the  sight  sec- 
tion alone.  As  test-section  power  increases  further 
(Photo  VI,  Fig.  7f),  the  voltage  trace  is  similar  to  that  of 
the  initial  description  of  the  STA  transition  except  that 
the  peak  voltages  are  not  full  deflections.  This  flow  is 
still  the  transitional  flow  from  slug  to  annular.   Physi- 
cally, it  appears  as  if  most  of  the  tube  consists  of  annular 
flow,  but  much  of  the  liquid  transport  is  still  in  the  very 
rapidly  moving  slugs . 

Annular  Flow  With  greater  power,  the  voltage  deflections 
become  smaller,  thus  indicating  that  the  annular  flow  region 
is  being  approached.  Annular  flow  was  defined  to  be  that 
flow  which  would  give  an  average  deflection  no  greater  than 
0.5  volts  (Fig.  7g)o  This  definition  did  allow  larger  de- 
flections to  occur,  but  the  frequency  of  these  large  deflec- 
tions was  greatly  reduced.   At  low  velocity  annular  flow 
was  well  defined  since  the  voltage  was  practically  zero 
(Photo  VII).   Photo  VIII  is  another  example  of  annular  flow. 
The  increased  velocity  of  the  vapor  core  caused  large  wave- 
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like  motions  on  the  surface  of  the  annular  film.  These 
waves  appeared  as  large  deflections  on  the  oscilloscope  . 
As  G  was  Increased,  these  waves  increased  in  frequency  and 
the  film  had  a  milky  appearance,  probably  due  to  the  in- 
creased entrainment  and  deposition  of  droplets.  Very 
little  could  have  been  said  about  the  interior  flow,  i.e. 
waves, if  only  photography  had  been  used  in  this  study. 
Wispy  Annular  Under  certain  conditions,  i.e.  T.   =  200°F, 
P  =  20-30  psia,  G  =1.6  -  2.0  x  106  lbm/hr-ft2,  a  varia- 
tion of  the  above  flow  regimes  was  noted.  As  power  was 
increased  beyond  bubbly  flow,  a  strong  bubbly  trace 
travelled  downward  from  (V,.)    to  0  volts.  Slug  flow  did 
not  appear  and  in  its  place  was  a  flow  which  can  be  called 
wispy  annular^'  or  froth*  '.   The  behavior  of  the  trace 
is  schematically  shown  in  Pigs.  7h,  7i,  and  7j  and  in  the 
corresponding  Photos  IX,  X,  and  XI.  As  the  void  fraction 
is  increased,  the  resistance  between  probe  and  wall 
becomes  greater  until  annular  flow  is  reached.  It  would 
be  possible  to  describe  many  flow  patterns  in  this  high- 
velocity  transition  between  bubbly  and  annular  flow. 
Froth  flow  would  probably  be  a  more  appropriate  descrip- 
tion of  the  beginning  of  this  transition,  while  wispy  annu- 
lar could  describe  the  latter  part  of  the  transition. 
Since  the  probe  indicated  that  the  transition  was  quite 
homogeneous  throughout,  it  appears  that  a  single  classi- 
fication is  adequate. 
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2.3.4  Charge  Carrying  Spray  Effect 

Photo  XII  again  shows  annular  flow,  but  the  oscillo- 
scope trace  is  radically  different  in  that  negative  traces 
up  to  one  volt  are  seen.  This  type  of  trace  was  not  ob- 
served in  high  pressure  studies^  '    '   and  initially  caused 
quite  a  bit  of  trouble  in  interpretation.  With  the  probe 
circuit  battery  removed  and  the  oscilloscope  connected 
directly  across  the  probe,  a  fluctuating  negative  voltage 
was  still  observed.  The  magnitude  and  frequency  of  this 
fluctuation  increased  as  the  flow  rate  was  increased,  even 
though  no  voltage  source  was  present.  A  possible  pickup 
effect  on  the  probe  from  the  voltage  across  the  test  tube 
was  discounted  because  the  generator  output  voltage  did 
not  have  such  large  fluctuations . 

The  explanation  of  this  negative  voltage  effect  is 
believed  to  be  as  follows.  As  annular  flow  is  approached, 
small  droplets  of  water  are  torn  away  from  the  liquid 
film  and  are  transported  as  a  spray  in  the  core.  By  some 
physical  effect  these  droplets  become  negatively  charged 
upon  detaching  from  the  bulk  flow,  thus  leaving  the  film 
positively  charged.  The  droplets  travel  down  the  tube  and 
rejoin  the  film  and  lose  their  charge,  escape  through  the 
exit  and  lose  their  charge  in  the  loop,  or  deposit  their 
charge  on  the  probe  .  This  last  possibility  appears  to  be 
the  source  of  the  negative  voltage  observed  on  the  scope. 
Professor  D.  P.  Keily^   '  has  observed  such  effects  with 
aerosol-type  flows  and  indicated  that  voltages  up  to 
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several  volts  could  be  observed  If  there  were  enough  drop- 
lets. The  mechanism  that  causes  the  detached  particle  to 
have  a  negative  charge  has  not  been  identified.  It  appears 
that  this  effect  is  quite  important  and  should  receive  more 
attention.  Vonnegut,  et  al.*  "'  indicate  that  such  charged 
water  droplets  are  produced  naturally  in  clouds  and  thunder- 
storms and  as  spray  from  lakes  and  oceans. 

This  negatively-charged  droplet  theory  is  substanti- 
ated by  several  observations.  When  the  water  was  particularly 
pure,  i.e.  3  megohm-cm  and  above ,  this  effect  was  quite 
strong.  However,  after  adding  salt  to  the  water  (30-40 
PPM),  the  negative  voltage  was  not  observed.  The  sketches 
k  and  1  in  Fig.  7  show  this  purity  effect. 

The  more  conductive  water  allows  more  of  the  negatively 
charged  particles  to  transfer  their  charge  to  the  liquid 
film  rather  than  depositing  it  on  the  probe .  This  also 
explains  why  the  Dynatech  workers,  who  used  a  similar 
probe,  did  not  observe  this  effect.  They  were  operating 
at  500  and  1000  psia  and  had  a  correspondingly  higher  temp- 
erature level.  The  conductivity  of  water  increases  greatly 
as  the  temperature  increases;  for  example,  the  conductivity 
of  pure  water  due  to  ionization  is  about  a  1000  times 
greater  at  546°F  than  at  77°F.^20'   Consequently,  this  con- 
ductive water  allowed  the  negative  charge  to  leak  into  the 
liquid  film  before  hitting  the  probe*  with  this  theory,  it 
would  be  expected  that  the  negative  voltage  would  increase 
with  increased  mass  flow  rate.  This  effect  was  also  obser- 
ved, as  stated  earlier. 
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Several  methods  were  attempted  to  eliminate  this  effect 
electronically  by  using  a  higher  battery  voltage,  and  by 
using  an  ac  high-frequency  voltage  source  .   The  first 
method  reduced  the  relative  voltage  of  the  effect  but  did 
not  eliminate  it.  By  means  of  filter  circuits  and  high 
frequency  ac,  it  was  hoped  to  eliminate  the  spurious  nega- 
tive charge  and  retain  only  the  desired  signals.  However, 
the  filters  and  bias  circuits  also  removed  the  desired 
fluctuations.  With  experience  it  was  found  that  proper 
interpretations  of  the  flows  could  be  made  even  though  the 
negative  voltage  was  present . 

2.3.5  Reproducibility  and  Reliability  of  Probe  Output 
About  a  dozen  runs  under  various  conditions  were 
repeated  for  a  second  time.  When  these  runs  were  compared 
to  the  original  ones,  it  was  seen  that  they  were  almost 
identical.   It  was  concluded  that  any  run  could  be  repeated 
with  the  expectancy  of  getting  very  similar  results. 

The  sight  section  observations  were  quite  helpful  in 
interpreting  the  probe  traces.  However,  when  the  bulk 
temperature  was  moderately  subcooled,  the  vapor  condensed 
rapidly  in  the  sight  section  and  it  was  necessary  to  look 
into  the  exit  of  the  heated  section  to  get  an  accurate 
visual  indication  of  the  flow  pattern.  When  the  void  frac- 
tion was  large,  the  flow  was  usually  so  disturbed  that  it 
was  difficult  to  visualize  the  flow  pattern.  In  both  of 
these  situations  the  probe  was  able  to  distinguish  the 
important  features  of  the  flow. 


2 .4  Test  Program 

After  the  probe  calibration  and  interpretation  runs 
were  completed,  the  actual  flow-regime  Investigations  were 
started.  The  following  table  lists  all  the  sets  of  runs 
completed  in  this  experimental  program. 

All  the  data  are  presented  on  flow  regime  maps.  The 
coordinates  of  mass  velocity  (G)  and  mass  quality  (x)  were 
chosen  because  it  appeared  that  they  were  best  suited  in 
presenting  the  various  effects  studied.   In  addition,  they 
are  usually  the  most  convenient  for  engineering  applica- 
tions. Researchers  in  adlabatic  two-phase  flow  frequently 
make  use  of  volume  flow  rate  and  volume  quality  coordinates. 
They  could  not  be  used  for  this  study  3ince  nonequllibrium 
effects  dominated  and  these  quantities  cannot  be  calculated. 
Several  runs  are  also  presented  in  terms  of  heat  flux  and 
quality  to  demonstrate  a  particular  effect. 

Due  to  the  manner  in  which  the  loop  operates,  it  was 
difficult  to  control  the  exit  pressure  very  accurately. 
Consequently,  pressure  comparisons  were  made  at  low 
pressures  where  the  exit  pressure  varied  from  20-35  psia 
and  at  high  pressures  between  90-100  psia.  These  small 
pressure  differences  did  not  affect  the  results  to  any 
appreciable  degree,  in  comparison  to  the  effects  due  to  the 
large  difference  between  high  and  low  pressures. 

2.5  Data  Reduction 

All  data  were  reduced  by  machine  computation.   The  IBM 
7094  Computer  at  the  MIT  Computation  Center  was  used.   The 
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program  is  described  in  Appendix  C.  The  heat  flux  was 
determined  by  the  formula 

The  equilibrium  exit  quality  is  given  by 
q/w  +  h,  -  h. 

*- tr1 -  (5) 


h 
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Chapter  III 
PRESENTATION  AND  DISCUSSION  OP  RESULTS 
3.1  Effects  of  Variables  on  Flow  Regimes 

Experiments  were  conducted  to  investigate  the  effects  of 
velocity,  inlet  temperature,  exit  pressure,  quality,  purity, 
length,  diameter  and  orientation  on  the  flow  regime  boundar- 
ies.  Figures  8  through  20  present  all  of  the  flow  regime 
data  in  terms  of  mass  velocity  versus  quality.   In  accordance 
with  the  discussion  in  Chapter  II,  the  following  flow  patterns 
were  identified:   bubbly,  bubbly-to-slug  transition,  slug, 
slug-to-annular  transition,  annular,  and  wispy  annular. 

It  is  seen  in  all  the  maps  that  the  bubbly  and  slug 
flows  occurred  primarily  in  the  subcooled  region,  while  full 
annular  flow  was  developed  at  2  to  4  percent  quality.  The 
two-phase  flow  is  therefore  in  a  highly  nonequilibrium  state 
and  a  considerable  volume  of  vapor  coexists  with  the  sub- 
cooled  liquid. 

From  all  the  maps,  it  is  observed  that  the  transitions 
occur  at  greater  subcooling  or  lower  quality  as  the  velocity 
of  the  flow  is  increased. 

The  influence  of  the  remaining  variables  on  the  flow 
regime  boundaries  can  best  be  illustrated  by  composite  plots. 
At  this  point  it  is  desirable  to  consider  only  the  major 
flow  regimes.   It  is  noted  that  the  transition  regions 
occasionally  cover  a  broad  enough  range  of  quality  so  that 
they  could  possibly  be  classified  as  flow  regimes.  However, 
practical  considerations  dictate  that  the  number  of  regimes 
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be  kept  to  a  minimum.   Boundary  lines  were  drawn  through 
the  center  of  the  transition  regions  and  these  lines  were 
then  transferred  to  the  appropriate  composite  plots.   These 
plots  are  presented  in  Pigs.  21  through  27. 

Figures  21  and  22  indicate  that  the  transitions  occur 
at  higher  quality  as  the  inlet  temperature  is  Increased 
from  100  to  200°F.   The  transition  bands  for  T,  =  200°F 
are  much  narrower  than  those  for  T.  =  100°F  (Figs.  10  and 
12) .   In  addition,  it  is  noted  that  slug  flow  was  not  ob- 
served  above  G  =  1.2  x  10u  lbm/hr-ft^  with  the  higher  inlet 
temperature.  These  latter  effects  were  also  observed  for 
the  smaller  diameter  test  section  (Fig.  19) ,  but  were  not 
observed  at  higher  pressure  (Fig.  14). 

At  T±  =  100°F,  at  30  and  100  psia,  particularly  at 
high  G,  the  transition  from  slug  to  annular  flow  was  notice- 
ably more  violent  causing  the  test  section,  the  connecting 
i*ubber  hoses,  and  piping  to  vibrate  considerably.  While 
testing  the  L/D  =  30  tube,  violent  vibrations  were  also 
observed  at  this  transition.  This  tube  burned  out  in  this 
region  in  spite  of  the  fact  that  the  heat  flux  was  not 
near  the  critical  flux  reported  by  Lopina.^   ' 

Pressure  had  a  considerable  effect  on  the  location  of 
the   flow  regime  boundaries.   Increasing  the  pressure  from 
30  psia  to  100  psia  shifted  the  flow  regimes  to  substan- 
tially lower  quality.   Figure  23  shows  this  pressure  effect 
at  Tx  =  100  and  200°F. 

In  order  to  eliminate  the  negative  voltage  problem 
previously  described,  salt  was  added  to  increase  the  conduc- 
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tivity  of  the  water.  The  effect  of  this  dissolved  Impurity 
on  the  flow  regimes  is  illustrated  in  Fig .  24.   The  BTS 
transition  for  very  pure  water  (2.5  megohm-cm)  occurs  at 
slightly  higher  quality  than  the  transition  for  the  1  and 
9  PPM  water.  There  was  no  discernible  effect  in  the  STA 
transition.  Conductivity  variations  between  1  and  3 
megohm-cm  produced  no  noticeable  shift  in  the  flow  regimes. 

Length-to-diameter  ratios  of  30,  60  and  90  were 
examined.  From  Fig.  25,  it  is  seen  that  the  flow  regimes 
are  quite  similar,  although  for  an  L/D  =  30  the  transitions 
are  at  slightly  lower  quality. 

For  the  same  T.  there  was  no  dlscernable  difference  in 
the  flow  regime  maps  for  the  two  diameters  (Fig.  26). 

In  Fig.  2?  it  is  seen  that  the  transition  lines  have 
the  same  general  shape  for  the  horizontal  and  vertical  test 

sections.  There  is  a  slight  shift  to  lower  quality  for  the 

6         2 
vertical  section.  At  velocities  below  0.4  x  10  lbm/hr-ft  , 

slight  stratification  was  visually  observed  for  the  bubbly 
and  BTS  transitions.  Since  the  vertical  test  section 
showed  only  a  minor  flow  regime  variation  from  the  horizon- 
tal section,  it  is  evident  that  this  stratification  did  not 
affect  the  probe  output  at  the  particular  position,  i.e. 
center,  that  the  probe  was  in.  A  modification  of  the 
probe  design,  whereby  the  probe  could  be  moved  radially 
across  the  tube  during  operation,  would  allow  examination 
of  the  flow  very  close  to  the  walls  and  thus  would  detect 
any  stratification. 
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3.2  Explanation  of  Results 

The  various  effects  noted  above  appear  to  be  very 
closely  related  to  incipient  boiling,  heat  flux,  and  tur- 
bulence of  the  flow.   To  explain  this  relationship,  the 
requirements  for  the  transition  to  slug  flow  are  first 

examined . 

(  21) 

Radovich  and  Moissisv   '  have  shown  that  the  frequency 

of  bubble  collisions  is  very  high  in  two-phase  bubbly  flow. 
A  relatively  small  number  of  these  collisions  results  in 
coalescence.  As  more  bubbles  coalesce,  the  bubbles 
increase  in  size  to  such  an  extent  that  slug  flow  is  obser- 
ved. The  time  that  it  takes  to  cause  the  first  identifiable 
existence  of  slug  flow  has  been  called  the  agglomeration 
time.  The  agglomeration  time  is  thus  determined  by  the 
collision  frequency  and  the  probability  of  coalescence  per 
collision.  The  latter  quantity  was  found  to  be  strongly 
dependent  on  the  purity  of  the  liquid.   Since  the  liquid 
used  in  this  study  was  maintained  at  relatively  constant 
purity  for  most  of  the  experiments,  the  agglomeration  time 
depended  primarily  on  the  collision  frequency.   A  decrease 
in  the  agglomeration  time  would  mean  a  shifting  of  the  BTS 
transition  to  lower  quality.  This  could  be  accomplished 
by  increasing  the  number  of  bubbles  and  by  causing  greater 
turbulence  in  the  flow. 

The  requirements  for  incipient  nucleation  in  a  forced- 
convection  surface-boiling  system  have  been  considered  in 

(22) 

detail  by  Bergles  and  Rohsenow.v   y  When  their  correlation 
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was  applied  to  the  present  conditions  (Appendix  D),  it  was 
found  that  incipient  boiling  occurred  at  greater  subcooling 
as  G  increased,  T1  decreased,  and  Pe  increased.  Further- 
more, it  can  be  seen  from  Figs .  28  to  31  that  the  heat  flux 
level  increases  as  G  increases,  T.  decreases,,  and  P 

increases.  As  the  heat  flux  increases  there  is  an  Increase 

(22) 

in  the  bubble  population.  Bowringv   '  has  concluded  that 

as  heat  flux  and  velocity  increase,  there  is  an  Increase 
in  the  subcooling  at  which  detatchment  of  the  bubbles  from 
the  heated  surface  occurs.   This  is  due  primarily  to  a 
general  increase  in  turbulent  mixing  under  these  conditions. 
This  enhanced  mixing  should  promote  the  agglomeration  of 
bubbles  . 

Thus,  on  the  basis  of  incipient  boiling,  heat  flux, 
and  mixing  considerations,  one  would  expect  a  shift  in  the 
BTS  transition  to  greater  subcooling  as  G  increases,  T. 

decreases,  and  P  increases.  These  trends  are  observed  in 

'      e 

the  present  data,  as  was  discussed  in  the  preceding  section. 
Since  the  heat  flux  level  increases  as  the  tube  length 
decreases,  the  shifting  of  the  BTS  transition  to  greater  sub- 
cooling as  length  is  decreased  is  also  accounted  for. 

Regarding  the  other  transitions,  it  is  reasonable  to 
expect  that  if  the  BTS  transition  occurs  at  lower  quality, 
the  STA  transition  should  also  occur  at  lower  quality.   In 
particular,  the  increased  bubble  agitation  at  the  higher 
heat  fluxes  would  be  expected  to  break  up  the  slugs,  thereby 
promoting  annular  flow. 
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The  limited  occurrance  of  wispy  annular  flow  can  be 
related  to  the  above  arguments.  It  Is  recalled  that  this 
regime  occurred  only  for  high  inlet  temperature  and  low 
pressure.  Under  these  conditions  the  heat  flux  level  is 
low  and  there  are  fewer  bubbles  and  also  less  tendency  for 
the  bubbles  to  agglomerate  into  slugs.  The  core  flow 
remains  bubbly  or  forthy  in  appearance  until  the  void 
fraction  is  so  high  that  the  spray  annular  configuration 
is  reached.   In  general  the  wispy  annular  regime  is  only 
observed  at  high  velocity  where  the  short  transit  times  are 
less  condusive  to  formation  of  slug  flow. 

At  this  point  it  is  worth  mentioning  that  Radovcich 
and  Moissis  predicted  that  the  BTS  transition  would  be  much 
more  rapid  for  pure  water.  As  discussed  in  the  preceding 
section,  the  present  tests  indicated  that  this  transition 
for  the  demineralized  water  occurred  at  about  the  same 
quality  as  that  for  the  water  with  dissolved  salt. 
Their  observation  was  for  adiabatic  flow;  thus,  with  heat 
addition  the  other  mechanisms  appear  to  overshadow  the 
impurity  effect . 

3.3  Interpretation  of  Critical  Heat  Flux  Data  by  Means  of 

Flow  Regimes 

The  present  flow  regime  observations  should  be  useful 
in  interpreting  heat  transfer  and  pressure  drop  data  for 
boiling  water  at  low  pressures.   As  an  example  of  the 
application  of  the  present  data,  one  of  the  important 
effects  in  bulk  boiling  burnout  will  be  considered. 
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Recent  studies  of  bulk  boiling  at  low  pressure  have 
confirmed  that  the  critical  condition  coincides  with  the 
disappearance  of  the  liquid  film  from  the  heated  surface.*   ' 
The  film  can  be  depleted  by  evaporation  and  entrainment,  or 
it  can  be  disrupted  by  bubble  nueleation.   In  addition  to 
these  local  effects  it  appears  that  unstable  upstream  flow 
patterns  can  produce  a  disrupt icn  of  the  annular  film 
further  downstream. 

Lopina^   '  found  that  the  burnout  heat  flux  increased 
as  the  inlet  temperature  increased.  Typical  data  illustra- 
ting this  effect  are  plotted  in  Fig.  33 o  The  flow  regimes 
for  these  conditions  are  superimposed  on  the  operating  lines. 

First  consider  operation  at  high  inlet  temperature  as 
shown  in  Fig.  33.  The  majority  of  the  tube  is  in  annular 
flow  with  only  moderate  increases  In  heat  flux.  The  slug 
flow  region  of  the  tube  becomes  increasingly  violent ,  but 
also  moves  toward  the  inlet,  as  the  heat  flux  is  increased. 
This  unstable  slug  flow  region  should  have  little  effect  on 
the  critical  heat  flux,  which  for  these  flow  conditions  is 
probably  caused  by  nuc lea tion° induced  breakup  of  the  film. 
Now,  as  the  inlet  temperature  is  decreased,  the  slug  flow 
region  remains  in  closer  proximity  to  the  exit  of  the  tube 
and  the  disturbance  created  within  this  region  would  be 
expected  to  propagate  downstream  and  affect  the  film  at  its 
thinnest  point,  which  is  at  the  exit  of  the  tube.  A 
significant  reduction  in  the  critical  heat  flux  is  thus 
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observed  as  the  Inlet  temperature  is  Increased. 

With  short  test  sections  the  operating  region  will  lie 
in  the  vicinity  of  zero  quality.   In  this  case  the  unstable 
slug  flow  regime  is  at  the  exit  to  the  test  section  and 
relatively  low  burnout  fluxes,  noted  in  the  present  L/D  = 
30  tests,  for  example,  would  be  expected. 
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Ghapter  IV 

CONCLUSIONS 
h  .1  Conclusions 

The  conclusions  of  this  investigation  can  be  summarized 
as  follows: 

1)  The  electrical  resistance  probe  is  exceptionally 
well  suited  for  flow  regime  studies .   It  is  simple,  inexpen- 
sive, and  provides  reliable  and  consistent  results  .  For 
low  pressure  studies,  where  nonequilibrium  effects  dominate, 
the  probe  is  more  reliable  than  an  exit  sight  section.. 

2)  The  interesting  phenomenon  of  the  spray  particles 
carrying  a  negative  charge  has  been  observed.  The  cause  of 
this  effect  is  not  yet  known,  but  various  researchers  are 
investigating  the  problem.  Although  it  influenced  the 
probe  output  under  certain  conditions,  it  did  not  inhibit  the 
flow  regime  analysis . 

3)  Plow  regime  observations  were  made  for  low  pressure 
water  over  a  wide  range  of  flow  conditions.  The  major  flow 
regimes  observed  were  bubbly,  slug,  annular,  and  wispy  annu- 
lar flow.   The  flow  regime  maps  show  that  bubbly,  transition 
to  slug,  slug,  and  transition  to  annular  flow  occur  in  the 
subcooled  region  under  nonequilibrium  conditions  while  a 
fully  developed  annular  flow  exists  by  the  time  2-4  percent 
quality  is  reached. 

4)  The  observed  effects  of  the  parameters  on  the 
flow  regimes  are  summarized  in  the  table  below. 
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Table  II 

Parameter 

Direction  of 
Change  in 
Parameter 

Direction  of  Shift  of 
Flow  Regime  Boundaries 

Mass  Velocity 

Increase 

towards  greater  subcooling 

Inlet  Temper- 
ature 

Increase 

towards  less  subcooling 

Pressure 

Increase 

towards  greater  subcooling 

Purity 

Decreased 

negligible  change 

Length 

Increased 

slight  change  towards  less  sub- 
cooling 

Diameter 

Decreased 

negligible  change 

Orientation 

Horizontal  to 

negligible  change 

vertical 

5)  The  influence  of  velocity,  temperature,  pressure  and 
length  on  the  flow  regime  boundaries  can  be  explained  in  terms 
of  incipient  boiling,  heat  flux,  and  turbulent  mixing  consid- 
erations. 

6)  The  present  observations  were  applied  to  bulk  boil- 
ing burnout  data.   It  was  found  that  the  effect  of  inlet  temp- 
erature on  burnout  can  be  related  to  unstable  upstream  flow 
patterns . 
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APPENDIX  A 

Burnout  Data 


Date:   19  and  21  July  1965 


Tl 

=  100°F 

l.d. 

=  0.242   in. 

L/D 

=  60 

Run 

Pe 
(psla) 

G  x  10 
(lbm/hr- 

-6 
ft2) 

X 

percent 

(q/A)cr  x  10-6 
(Btu/hr-ft2) 

1 

27.0 

0.5 

33.0 

0.962 

2 

28.0 

0.4 

37.0 

0.843 

3 

43.5 

0.8 

30.6 

1.51 

4 

34.0 

0.8 

19.0 

1.12 

5 

37.5 

1.2 

11.6 

1.36 

6 

36.0 

1.2 

10.9 

1.32 

7 

51.0 

1.6 

IO.76 

1.81 

8 

45.O 

1.6 

8.35 

1.64 

9 

48.0 

2.0 

4.15 

1.71 
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APPENDIX  B 
Photographic  Procedures 

For  the  photographic  study,  simultaneous  pictures  were 
taken  of  the  oscilloscope  trace  and  of  the  exit  flow  in  the 
glass  sight  section.  This  was  accomplished  by  shooting  both 
pictures  manually  at  the  same  time  at  the  desired  flow 
conditions . 

The  Polaroid  Oscilloscope  Camera  for  the  Hewlett-Packard 
Company  Model  130B/BR  Oscilloscope  was  used.  The  film  was 
Polaroid  ASA  3000  film  and  the  camera  settings  were  f  2,8 
and  speed  of  1/30  sec . 

The  sight  section  was  normally  illuminated  by  a 
General  Radio  Co  ,  Strobotac  .  The  Strobotac  was  always  off 
when  actually  taking  photos. 

For  the  sight  section  photos,  Polaroid  ASA  200  Film  was 
used  in  order  to  obtain  better  quality  pictures.  A  General 
Radio  Co.  Type  1530-A  Microflash  was  used  to  illuminate  the 
exit  test  section.  The  2-microsecond  light  pulse  effectively 
stopped  the  motion  of  the  flow  in  the  sight  section.  Various 
backgrounds  for  the  sight  section  were  tried  with  black 
velvet  giving  the  best  results .  The  Polaroid  camera  was 
mounted  on  a  tripod  and  the  +2  and  +4  Polaroid  close-up 
lenses  were  used.   The  camera  settings  were  f.8  and  1/30  sec. 

Below  is  a  sketch  of  the  lighting  setup  viewed  from 
above.  The  microflash  was  approximately  12  inches  below  the 
test  section  and  was  inclined  so  that  the  light  was  pointed 
directly  at  the  sight  section. 
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APPENDIX  C 
Computer  Program  for  Data  Reduction 

The  experimental  data  obtained  was  reduced  in  accor- 
dance with  the  equations  in  Section  2,5  by  a  program  written 
in  FORTRAN  and  run  on  the  MIT  Computation  Center  IBM  7094 
Computer. 

The  input  data  consisted  of  the  run  identification 
number,  the  exit  pressure,  voltage  across  the  test  section, 
shunt  voltage  (current  through  test  section),  enthalpy  of 
the  incoming  fluid, and  the  mass  flow  rate. 

The  program  calculated  the  rate  of  heat  transfer,  heat 
flux,  rate  of  heat  transfer  per  pound  mass  of  water,  enthal- 
py of  saturated  liquid,  latent  heat,  and  quality. 

Below  is  a  list  of  the  variables  as  used  in  the  program 

and  their  physical  counterpart . 

PROGRAM  VARIABLE  PHYSICAL  COUNTERPART 

I,  J,  K  Numeral  run  identification  numbers 

P  EXIT  P  -  the  pressure  at  the  test  section 

exit  (psia) 

E  V  -  voltage  across  test  section  (volts) 

CUR  I  -  shunt  voltage  (mv) 

HIN  h.   -  enthalpy  of  incoming  subcooled 

in  water 

W  w  -  mass  flow  rate  (lb/hr) 

A  A  -  heat  transfer  area  (Changed  this 

value  for  each  L/D  and  D)  (ft2) 

HADDED  q  -  rate  of  heat  flow  (Btu/hr) 

o 
HTPLUX  q/A  -  heat  flux  ( Btu/hr -ft  ) 
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BTULBM 
ENFPX 


q/w  -  rate  of  heat  flow  per  lbm  of 
water  (Btu/hr-lbm) 

hfe  -  enthalpy  of  saturated  liquid  at 
exit  (Btu/lbm) 


h 


fg  -  latent  heat  at  exit  (Btu/lbm) 


ENFGPX 

X  x  -  quality 

The  enthalpy  of  saturated  liquid  (hf)  as  a  function 

(  25) 

^'   by  a   power   series 

For   P  <   200   psia 


of   pressure  was   obtained  by  Todreas 
fit   of  the   Keenan  and  Keyes'    values 


hf(P)   =    Z  a   P1 

i=o  x 


where 


aQ  =  1. 1222734(10^) 


a±   =  6.3204790(10u) 


-1 


a      =  2.1239014(10"7) 

a6  =  -1.0067598(10~9) 

a?   =  2.9480958(10"12) 


a8  =  -4,8430355(10 


-15 


a2  =  -1.4742752(10~J-) 

a3  =  2.5403593(10~3) 

a^  =  -2.8788220(10"5) 

The  latent  heat  (h-  )  was  obtained  by  fitting  the 
Keenan  and  Keyes'  values  on  a  power  series.  For  P  <400  psia, 


3.4075972(10"18) 


h  (P)  =  Z  a  P1 
1&  i=o  1 


•1.2499843(10"9) 


12- 


aQ  =  9-9704457(10^) 
a±   =  -2.1762821(10°) 
a2  =  1.9558791(10"2) 
a3  =  -1.2610018(10"^) 
a^  =  5.0308804(10"7) 

The  values  of  the  heat  transfer  area  for  the  differ- 
ent tubes  are: 


a6  =  1.9322768(10"-^) 
a   =  -1.8023251(10"15) 
ag  =  9.2680092(10"19) 
a9  =  -2.0148643(10"22) 
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TUBE 

A  (ft2) 

0.242  in.  i .d .  L/D 

=  30 

0.038274 

0.242 

60 

0.076548 

0.242 

90 

0.114822 

0.094 

60 

0.0154 

The  output  data  consisted  of  all  the  input  variables 
and  the  computed  quantities. 

The  program  used  for  the  L/D  =  60,  D  =  0 ,094-in.  test 
tube  computation  follows. 

COMPUTATION  PROGRAM 

1  READ  2,  I,  J,  K,  PEXIT,  E,  CUR,  HIN,  W 
A      =  .0154 

HADDED  =  205.36  *  E  *  CUR 
HTFLUX  =  HADDED  /  A 
BTULBM  =  HADDED  /  W 
P      =  PEXIT 

5  ENFPX  =  1.1222734E2  +  6.3204790*P  -  1 .4742752E-1*(P**2) 

1  +2.5403593E-3  *(P**3)  -2 .8788220E-5*(P**4)  +2.12390l4E-7*(P**5) 

2  -1.0067598E-9*(P**6)  +2 .9480958E-12*(P**7) 

3  -4.8430355E-15*(P**o)  +3 .4075972E-l8*(P**9) 

6  ENFGPX=9 .9704457E2  -2.1762821*P  +1 .955o791E-2*(P**2) 

1  -1.26l00l8E-4*(P**3)  +5.0308804E-7*(P**4)  -1 .2499843E-9*(P**5) 

2  +1.9322768E-12*(P**6)  -1 ,8023251E-15*(P**7) 

3  +9.2680092E-19*(P**8)  -2  ,0l48643E-22*(P**9) 
X  =  (HIN  +  BTULBM  -  ENFPX )/ENFGPX 

PRINT  3,  I,  J,  K,  PEXIT,  E, 
1   CUR,  HIN,  W,  HADDED,  HTFLUX,  BTULBM,  X,  ENFPX,  ENFGPX 
GO  TO  1 

2  FORMAT  (315,5E10.2) 

3  FORMAT  (7H0RUN  NOI3,  1H-,  13,  1H-,13,5X, 

1  6HPEXIT=E15.6,5X,2HE=E15.6,5X,4HCUR=E15.6/5X, 

2  4HHIN=E15.6,5X2HW=E15.6,5X,7 

3  EIHADDED=E15  •  6/5X, 7HHTFLUX=E15  .6, 5X, 

4  7HBTULBM=E15  •  6 ,  5X,  21IX=E15  . 6/5X 

5  6HENFPX=E15.6,5X,7HENFGPX=E15 .6) 
END 
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APPENDIX  D 
Incipient  Boiling  Analysis 
The  effects  of  velocity,  temperature,  and  pressure  on 
the  Incipient  boiling  heat  flux,  (q/A).,  can  be  examined  by 
means  of  the  correlations  given  in  Ref .  22.  This  correla- 
tion for  water  is: 

2,30 


(q/A)1    =  15.60    P1J56(Tm-t/  (Dl) 

For  a  uniformly  heated  tube  the  heat  flux  can  be  calculated 
from: 

(q/A)  =  h(Tw-Tb)  =  h  (Tw-Ts)  +  (T^)  (D2) 

The  saturation  temperature  is  known  from  the  exit  pressure. 
The  bulk  temperature  is  found  by  the  First  Laws 

*   =  °p<Tout-Tin)w  <D3) 

which  leads  to: 

Tout  =  (^A)^  +  Tln  =  Tb  •  W 

Since  the  probe  is  at  the  exit,  only  the  conditions  at  the 

exit  are  of  interest  and  therefore  TK  =  T  .  .  The  heat 

b  out 

transfer  coefficient,   h,    is   found  using  McAdams1    correla- 
tion:<26> 

(f)   =0.023   (fr8(^K4  (D5) 

where  the  properties  are  evaluated  at  the  bulk  fluid  tempera- 
ture .  For  the  actual  numerical  analysis,  the  values  of  these 

(  07) 

properties   came   from  ElWakil.v    ' '     By  assuming  a   value   of 
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(q/A),  T,  and  hence  h  can  be  calculated.   Then  using  (D2), 
(T  -T  )  is  found.  The  values  of  (q/A)  and  (T  -T  ),  for 

W   S  W   s 

various,  G,  T.,  and  P  can  then  be  found  and  plotted. 
Equation  (Dl)  for  the  various  pressures  is  then  superimposed 
on  these  plots.  The  intersections  of  the  (q/A)  vs.  (T  -T  ) 
curve  and  the  (q/A)1  curve  is  the  point  at  which  Incipient 
boiling  is  predicted. 

This  analysis  was  carried  out  for  the  conditions  that 
were  presented  in  Figs.  28  through  31,  i.e. 


G  x  10 
lbm/hr- 

-6 
ft2 

°F 

Pe 
(psia) 

a) 

0.4 

100,  150,  200 

30 

b) 

0.4 

100,      200 

100 

c) 

2.0 

100,  150,  200 

30 

d) 

2.0 

100,      200 

100 

In  Fig.  32a  is  a  representative  plot  for  the  conditions  in 
(c).  To  visually  relate  the  predicted  (q/A).  to  (x),  opera- 
ting lines  obtained  from  the  formula: 

q/w  +  h,   -  hr 
x  =  hln    f  (D6) 

nfg 

are  plotted  on  coordinates  of  (q/A)  vs  x.  The  (q/A).  found 
from  the  analysis  are  plotted  on  these  operating  lines  and  are 
connected  by  a  smooth  line.  This  line  then  is  the  locus  of 
predicted  (q/A).  for  various  inlet  temperatures,  This  plot 
for  conditions  (c)  above  is  presented  in  Fig.  32b.  Examin- 
ing Fig.  32b  it  is  seen  that  incipient  boiling  occurs  at 
greater  subcoolings  for  lower  T . .  From  the  analysis  of 
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conditions  a,  b,  d,  it  was  also  found  that  incipient  boiling 
occurs  at  greater  subcooling  for  larger  pressures  and,  keeping 
P  and  T.  constant,  at  higher  G. 
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FIG.    2      DETAILED  VIEW   OF   TEST   SECTIONS 
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FIG.    3a      PROBE   CIRCUIT   DIAGRAM 
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FIG.  3b   CHARACTERISTICS  OF  PROBE  CIRCUIT 
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PHOTO  II 


P~30,  G=2.0xl0  ,  T±=200,  x=-2.8lSi 
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PHOTO   III 
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PHOTO   IV 


P~30,    0=0.4x10    ,    T±=150,    x=-.20l£  P~30,    G=2.0xlOb,    T±=100,    x=-1.43£ 

SLUG  SLUG 

FIG.    4      COMPARISON    OF    PROBE   AND   SIGHT-SECT  SERVATIONS 

(BUBBLY,    BUBBLY -TO -SLUG,    SLUG) 
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P~30,    G=0.4xl06,    T.=200,    X=H   • 

P~30,    G=1.2xl0    ,    T±=_       . 

ANNULAR 

ANNULAR 

FIG.  5   COMPARISON  OF  PROBE  AND  SIGHT-SECTION  OBSERVATIONS 
(SLUG -TO -ANNULAR,  ANNULAR) 


-52- 


1.0 


IT' 


PHOTO  IX 


;mvk 


zee 


,b 


0.04 


P-v30,  G=2.0xl0  ,  T±=200,  x=-0.T6c,b 
WISPY  ANNULAR 


PHOTO  X 

P~30,  G =2. 0x10°,  T.=200,  x=-0.66^ 


WISPY  ANNULAR 
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PHOTO  XI 
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P~30,  G=2.0xl0  ,  T.=200,  x=4.69$ 
WISPY  ANNULAR 


PHOTO  XII 

P~30,  G=0.6xl0°,  Ti=200,  x=ll. 

ANNULAR  WITH  NEGATIVE  CHARGE  EFFECT 


FIG.  6   COMPARISON  OF  PROBE  AND  SIGHT -SECTION  OBSERVATIONS 
(WISPY  ANNULAR) 
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